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Considerable interests on the charmonium states in the medium [1] have been triggered ever 
since the proposal of considering the dissolution of the charmonium states as a signal of the forma- 
tion of QGP [2]. However, the understanding of the fate of the charmonium states in the medium 
still remains very difficult. The potential model approach at finite temperature is still under scrutiny 
since the potential used in the Schrodinger equation is not rigorously defined [3]. Remarkable 
progress has been made recently in the heavy quark effective theory at finite temperature [4], where 
the potential can be rigorously defined. However, it requires scales to be hierarchically ordered. 
More recently, a new approach based on the path integral formalism has been proposed [5]. It 
needs further research since the effects of the medium on the heavy quark in this approach is mod- 
eled with only Coulomb interaction. First principle calculations in lattice QCD are thus crucially 
needed to determine the properties of the charmonium states in the hot medium. 

With the lattice QCD approach [6], the properties of the charmonium, which can be directly 
seen from the spectral function, are contained in the Euclidean time correlation functions. The 
extraction of spectral functions from correlators is rather difficult due to the limited number of 
points in temporal direction required to perform an analytic continuation from imaginary to real 
time. In this work, we present the analysis results of charmonium properties at both the correlator 
and spectral function level. Previous work has been reported in Ref. [7, 8]. 

The charmonium correlators at vanishing momentum are calculated using: 

G h (t,T) = j £(J h (t,x) 4(0,0) ) T . (1) 

x 

Jh is a suitable mesonic operator, here we consider a local operator of q(z,x)THq{z,x), where 
r# = JuYs for vector (V,,) and pseudo-scalar (PS) channels, respectively. The temperature T is 
related to Euclidean temporal extent aN T by T = 1/ (aN x ), where a is the lattice spacing. Through 
analytic calculation, the correlation function can be related to the spectral function as the following: 



p co 

G H {z,T) = dcoo H ((0,T)K(z,T,co), 
Jo 



(2) 



where the kernel K is given by K(z,T,(o) = cosh(eo(T — ^))/sinh(^). The spectral function 
o(co) contains all the information of the hadron properties in the medium and is the key quantity to 
be investigated. For instance, the dissociation temperature can be read from the deformation of the 
spectral function and the heavy quark diffusion constant D relates to the vector spectral function as 

3*00 03 

where Xoo is the quark number susceptibility. 

Inverting Eq. @ to obtain the spectral function at finite temperature is hampered mainly by 
two issues: the temporal extent is always restricted by the temperature, az < l/T; the spectral 
functions we want to have should be continuous and have a degree of freedom of ^(1000) but 
the correlators are calculated in the discretized time slices with limited numbers, typically £?(10), 
which makes the inversion ill-posed. Thus the normal % 2 fitting would be inconclusive. Maximum 
Entropy Method (MEM) is currently one of the best tools in the literature to solve the problem [9]. 
It is based on the Bayesian algorithm and requires the prior knowledge of spectral functions as an 
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J8 a [fin] a _1 [GeV] L a [fin] c sw * A^x/V T T/T c iV, 



conf 



7.793 0.010 18.974 1.33 1.310381 0.13200 128 3 x 96 0.73 234 

128 3 x 48 1.46 461 



128 3 x 32 2.20 105 
128 3 x 24 2.93 81 



Table 1: Lattice parameters. 



input. The only parameter in the MEM analysis is the so called "default model", which provides 
the prior information of the spectral function, e.g. the spectral function should be positive-definite. 

The standard Wilson plaquette action for the gauge field and the non-perturbatively G(a) im- 
proved clover fermion action for charm quarks are implemented in the simulation. The mass of 
vector meson is tuned to the physical J/y mass. We measured correlation functions on very fine 
quenched lattices with a relatively large size of 128 3 x 96, 128 3 x 48, 128 3 x 32 and 128 3 x 24 at 
0.73 T c , 1.46 T c , 2.20 T c and 2.93 T c , respectively. The lattice parameters are shown in Table RL 

We first analyze the temperature dependence of charmonium states at the correlator level. One 
constructs [10] 

r , , i\ cosh (o>(t- 1/27/)) 

G rec (z,T;T' )= / dOCT fl>,7") ) ) 1 >\ (4) 

Jo sinh(ft)/2r) 

to study the difference of the spectral function at temperature T and T'. The deviation of G(z,T) 
from G Tec (z,T;T') indicates any modifications of the spectral function at temperatures T from the 
one at temperature T'. One normally needs a technique to obtain the spectral function a (ft), 7/') at 
a reference temperature T' and consequently the evaluation of Eq. (Q) suffers from the uncertainty 
of the determination of the spectral function brought by the certain technique. We found a useful 
exact relation, which is a generalization of Ref. [1 1], as follows: 

cosh[a)(T-/vy2)] = cosh[fl)(T'-/V;/2)] 

sinh(fl)/V r /2) ~ T , =T; V+=^ sinh(fl)/V£/2) ' M 

where T' = {aN' x )- x , T = (aN % )-\ x' 6 [0, A^- 1], T€[0,/V T -1], N' r = mN x , meZ+./V T and 
N' T are the number of time slices in the temporal directions at temperature T and T', respectively. 
T denotes the time slice of the correlation function at temperature T while z' denotes the time slice 
of the correlation function at temperature T . The sum of z' on the right hand side of Eq. (||) starts 
from z' = z with a step length N T to the upper limit N' Z — N T + z. After putting a(ft), T') into both 
sides of the above relation and performing the integration over ft), one immediately arrives at: 

N' T -N T +T 

G Tec (z,T;T')= £ G(z',T'), (6) 

t'=t; x'+=N x 

which shows the evaluation of G iec (z,T;T r ) can be done directly from the correlator G(z',T') at 
T 1 . In what follows, we suppress T' in G iec (z, T; T'). 

The ratios G(z,T)/G Tec (z,T) are shown in Fig. [IJ The G rec are obtained from the correlator 
data at 0.73 T c through Eq. @. In the left plot of Fig. [IJone can see that the ratios in the PS channel 
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Figure 1: The ratio G(t,T)/G kc {t,T) in PS (left) and Vu (right) channels at T = 1.46, 2.20 and 2.93 T c . 



at all temperatures decrease monotonically with increasing distance. The temperature effects set in 
at a shorter distance at a higher temperature. At the largest distances, the ratios deviate from unity 
at around 5%, 8% and 12% at 1.46, 2.20 and 2.93 T c , respectively. The results for the Vu channel 
are shown in the right plot of Fig. [j]. Different from the case in the PS channel, the ratios increase 
monotonically with increasing distance. At the largest distances, the ratios have a larger deviation 
from unity compared with the results in the PS channel. The deviation at the largest distance at 
1.46 T c is around 16% and it is comparable to the values at higher temperatures. 
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Figure 2: The ratio G sub (T,r)/G r s e u c b (T,r) (G m (z,T)/G^ (t,T)) in PS (left) and Vu (right) channels at 
T = 1 .46, 2.20 and 2.93 T c . The "diff" and "sub" stand for the results of the ratio G diff /G™ and G sub /G™ b . 

The large temperature dependence of G(r,T)/G rec (r,T) in the Vu channel could be due to the 
possible zero mode contributions. To suppress the zero mode contribution we evaluate the ratio of 
the differences of the neighboring correlators to the difference of the corresponding reconstructed 
correlators [12] 



G dlft (T,r) _ G(z,T)-G(z+l,T) 



(J) 



G^(r, T) G rec (T, T) - G rec (T+ 1, T) ' 
which equals the ratio of the time derivative of the correlators to the time derivative of the recon 
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structed correlators at X + 1 /2. One can also check the ratio of midpoint subtracted correlators [13] 

G™ b (T,T) G(z,T)-G(N t /2,T) 
Gg*(T, T) G vec (T, T) - G rec (jV T /2, T)' 

In those two ratios a % independent constant cancels in the correlator. 

The results for G sub /G™ b and G dm /G?^ in PS (left) and V ti (right) channels are shown in 
Fig. |[ The ratios G sub /G™c and G dlft /G d gf give similar results at all the distances. Seen from the 
right panel of Fig. ||the magnitude of the measured vector correlator to the reconstructed correlator 
reduces dramatically after the implementation of the difference of neighboring correlators (Eq. ([7])) 
and mid-point subtracted correlators (Eq. (||)). At 1.46 T c the ratio is more or less unity at all 
distances, at 2.20 T c and 2.93 T c the ratio becomes even smaller than unity at large distances. The 
zero mode contribution accounts for most temperature dependence of the ratio G / G rec at least at 
1 .46 T c seen in Fig. [j]. The deviations of the ratios from unity in the PS channel shown in the right 
panel of Fig. ^| are also reduced. However, the effect is not as strong as that in the Vu channel and 
the values at the largest distance are shifted up only about 3% at both 1.46 T c and 2.20 T c and about 
6% at 2.93 T c compared with the results in Fig. [|. Comparing the results for the Vu channel with 
those for the PS one in Fig. |, we find the ratios G sub /G™ b (G diff /G d if ) in these two channels have 
similar behavior at all distances at two higher temperatures 2.20 T c and 2.93 T c . However, they 
differ at 1.46 T c . The phenomenon we observe here might suggest that J/y could survive up to 
1.46 T c and starts to melt at 2.20 T c , and tj c might be melted already at 1.46 T c . 

One has to note that the comparison of the measured correlator with the reconstructed correla- 
tor can only give a rough idea of the magnitude of any medium effects at a certain temperature. To 
really explore the properties of the charmonium states, one has to go to the spectral function level. 
Thus it is crucial to extract the spectral function from the correlators using MEM. 

In the MEM analysis, we use number of points in the investigated energy region N® = 8000, 
the minimum energy a(O m i„ = 0.000001 and implement the modified kernel K* = tanh(ft)/2) K [8, 
14] to explore the low energy behavior of spectral function [15]. The default model we used in the 
PS channel is a normalized free lattice spectral function and in the Vu channel is a normalized free 
lattice spectral function with a transport peak modeled by a Breit-Wigner distribution at small co. 
Here we show the spectral functions in PS and Vu channels with statistical uncertainties 1 . 

The statistical error is obtained from the Jackknife method. Quite often the statistical error is 
given on the integral of the spectral function a (a) over a certain CO region in the spectral function 
plot [10, 17]. However, it is not straightforward to get a feeling of how big the error is on the 
spectral function itself. Here we rather calculate the Jackknife error on each point of the spectral 
function. We show the results for PS (left) and Vu (right) channels in the intermediate frequency 
region in Fig. |3] and show the transport part of spectral functions in the V„ channel in Fig. f|. In 
both figures, shaded areas denote the Jackknife errors and solid lines inside the shaded areas are 
the mean values of spectral functions. 

From the left plot in Fig. || one can see that at 0.73 T c the spectral function in the PS channel has 
large uncertainties in the amplitude at the point which corresponds to the ground state peak location 
in the mean spectral function. However, even at the lower end of the error bar, the amplitude is still 

'The systematic error analyses have been performed in Ref. [16], which do not change the general results shown in 
Fig. | and Fig. Q 
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Figure 3: The statistical errors for the output spectral functions in PS (left) and V„ (right) channels at all 
available temperatures. The shaded areas are the errors of the output spectral functions. The mean values 
are the solid lines inside the shaded areas. 




0.75 1 1.25 1.5 1.75 2 2.25 2.5 2.75 3 



Figure 4: Left: the statistical errors for the transport part of output spectral functions in the V„ channel at all 
available temperatures. The shaded areas are the errors of the output spectral function. The mean values are 
the solid lines inside the shaded areas. Right: the quark number susceptibility Xoo/T 2 versus T /T c . 



larger than the peak amplitudes at the higher temperatures within the errors. The peak location of 
the ground state peak at 0.73 T c might be shifted to a lower energy of CO « 3 GeV or to a higher 
energy at CO ~ 3.6 GeV. In the latter case, the peak location would have the same peak location as 
the spectral function at 1.46 T c but with a much larger amplitude and smaller width. At 2.23 T c 
there is hardly a peak structure within the statistical errors. At 2.93 T c the spectral function flattens. 
Thus this picture suggests r\ c is "partly" melted at 1.46 T c and dissolves at higher temperatures. In 
the right plot of Fig. [3], we focus on the resonance part of the spectral function in the Vn channel. 
One sees that the peak location of the spectral function at 1 .46 T c does not have an overlap with the 
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peak location of the spectral function at 0.73 T c . The amplitudes between these two differ a lot. At 
both 2.20 T c and 2.93 T c there are hardly any peak structures and at 2.93 T c the spectral function is 
flattened. This picture indicates 7/i/^is already melted at 1.46 T c . 

The statistical uncertainties of the transport peaks in the V,, channel are shown in the left 
plot of Fig. Q The amplitude of the transport peak at (O = gives the value of the heavy quark 
diffusion constant. The uncertainties of both, amplitudes and widths of the peak, are relatively 
small. Through Eq. <j3j> and Xoo/T 2 shown in the right plot of Fig. f|, we get DT very roughly to be 
0.28 at 1.46 T c and find DT increases with increasing temperatures. The precise determination of 
DT needs further detailed study. 

In summary, our analysis suggests that J j\\f is melted already at 1.46 T c and r\ c starts to dis- 
solve at 1 .46 T c and does not exist at higher temperatures. We also identify the transport contribu- 
tion at the spectral function level for the first time. 
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